The A(1232) axial charge and form factors from lattice QCD 
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We present the first calculation on the A axial-vector and pseudoscalar form factors using lat- 
tice QCD. Two Goldberger-Treiman relations are derived and examined. A combined chiral fit is 
performed to the nucleon axial charge, N to A axial transition coupling constant and A axial charge. 



Introduction. Numerical solution of Quantum Chro- 
modynamics (QCD), the underlying theory of the strong 
interactions, has proved a very successful approach in 
providing a theoretical understanding of baryon struc- 
ture. During the last few years, simulations of the dis- 
cretized theory known as lattice QCD have included dy- 
namical quarks with near to physical mass values [ij. A 
success of these recent simulations is the calculation of 
the low-lying hadron spectrum ^ iJ] showing agreement 
with the experimental values. The lattice set-up can also 
be applied to compute quantities that are not known ex- 
perimentally. In this work we report the first calcula- 
tion of the A axial- vector and pseudoscalar form factors 
(FFs). 

Understanding the structure of the A resonance has 
great relevance to nuclear phenomenology. The A is a 
rather broad resonance close to the ttN threshold. It 
therefore couples strongly to nucleons and pions mak- 
ing it an important ingredient in chiral expansions . 
The A baryon resists experimental probing due to its 
short lifetime {^^ 10~^^ s) 0, [l^l- Its axial charge and 
TT-A coupling constants that are needed as input in chiral 
Lagrangians are difficult to measure. Baryon chiral ex- 
pansion calculations that include the A explicitly follow 
one of two strategies as far as the determination of these 
parameters is concerned. The first is to relegate the ax- 
ial charge to one of many fitparameters, and fit using 
lattice [5|, |ll| , experimental (8| , or partial- wave calcula- 
tion data [7|. The second is to use estimates based on 
phenomenology such as the relation between the nucleon 
axial charge qa that is well measured and the A axial 
charge, which can be derived from the large- iVc hmit [l2j 
or SU (4) symmetry [l^ . The Goldberger-Treiman rela- 
tion is then used to get the effective ttAA coupling. 

Quite recently, groups have calculated A axial charge 
through QCD sum rules [lj| and xPT 15|, and both have 



noted the lack of an explicit lattice calculation of this 
quantity. First-principles lattice QCD calculations can 
probe the structure of the A and indeed recent studies 



have produced calculations of the ttA^A coupling [16l |17| 
and the electromagnetic form- factors of the A [18| . Using 
our lattice QCD formulation, we are then well-positioned 
to calculate the axial charge of the and the effective pion- 
A coupling, GttAA, as well as examine the Goldberger- 
Treiman relations as a way to relate the A axial charge 

to GttAA- 

In this letter we present the first lattice calculation of 
the axial- vector and pseudoscalar form factors of the A 
baryon. At non-zero momentum transfer , we find a 
second pseudoscalar form-factor, yielding a second effec- 
tive coupling constant, iJTrAA- These calculations lead 
to two Goldberger-Treiman relations, which are indeed 
satisfied by the lattice results. 

Axial-vector and Pseudoscalar Matrix Element. We 
consider the matrix element of a current X between A+ 
states 

(A(p/,S/)|X|A(p,,Sj)) =U„{pf,Sf) [O^Y^ Ur{p^,Si), 

where u„ denotes the Rarita-Schwinger vector-spinor 
(cr, T are Lorentz indices), and pi and pf are the initial 
and final momenta of the A. For the axial-vector cur- 
rent A^j^ix) ~ 'ijj{x)^f^"f^^'ip{x) the matrix element can 
be written in terms of four Lorentz-invariant FFs, labeled 
gi, gs, hi and /13: 
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where q = Pf — Pi is the momentum transfer. For the 
pseudoscalar current P°'[x) = 'ip{x)^z^ip{x), the matrix 
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element can be written in terms of two FFs, to be defined 
below in relation to the Goldberger-Treiman relations. 

Lattice Evaluation. Axial form factors on the lattice 
are extracted in a standard way from the three-point 
function 

(Gt^'^(</,i;p/,p.;rp))= e-w-^e+'''-^ r^" 

X2. Xi 

{n\T [xr(x2,i2)4(xi,ii)xl^(o,o)] \n){2) 

where XA is an interpolating operator creating a state 
with the quantum numbers of the A"*" [l^ and Tp is a set 
of projectors, given by r4 = i(l-|-74) and = ir4757fe. 
A similar three-point function with is required 

for the extraction of the pseudoscalar FFs. Technically 
these are evaluated via the sequential inversion through 
the sink [18| at fixed sink time-slice tf, while the 



V 



stat. m,r (Gov) ruN (GeV) niA (GeV) 



and operator insertion is supplied at all intermedi- 
ate t-slices {0 < t < tf) and Fourier-transformed for all 
momenta q at a small extra CPU cost. 

The kinematics arc fixed to a static A sink (py = 
0,q = — Pi). Denoting for convenience Monte-Carlo av- 
erages G^,(r„q,i) = (G^;^^(i;, t ; 0, -q; Fp)) for X = 
or P^, we construct the optimal ratio of three-point 
to two- point functions 

/ GA(r4,-q,tj-t)GA(r4,o,t)GA(r4,o,t/) 
V GA(r4,o,t/-i)GA(r4,-q,<)GA(r4,-q,</) 



n;);(rp,q) 
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where GA(r^,p,t) is the A propagator of momentum 
p. This ratio eliminates unknown field renormalization 
constants and leading time dependences and tends to a 
constant at large Euclidean time separations tf — t and 
t. A careful optimization in the space of the source-sink 
Lorentz indices a, r, p is required and only two linear com- 
binations of sequential sources suffice to provide all four 
axial and two pseudoscalar FFs. 

Smearing techniques are implemented resulting in sat- 
isfactory suppression of excited state effects allowing the 
source-sink distance to be fixed at about 1 fm 18|. Lat- 
tice computations of the matrix elements of the axial- 
vector and pseudoscalar currents for all transition mo- 
menta vectors q contributing to a given value of = 
— (pf — PiY are simultaneously analyzed and the over- 
constrained system determines the form factors through 
a global minimization. We note that 0(500) lattice 
measurements are involved in the extraction of the form 
factors for Q^-values up to ~ 2 GeV^. 

The parameters of the lattice ensembles used in this 
calculation are given in Table H] The quenched Wilson 
fcrmions gauge configurations enable the extraction of 



Quenched Wilson fermions 
P = 6.0, = 2.14(6) GeV 



32'' X 64 200 
32^ X 64 200 
32^ X 64 200 



0.563(4) 
0.490(4) 
0.411(4) 



1.267(11) 
1.190(13) 
1.109(13) 



1.470(15) 
1.425(16) 
1.382(19) 



Mixed action, a"^ = 1.58(3) GeV 
Asqtad (am„,d/= = 0.02/0.05), DWF (am„,d = 0.0313) 
20^ X 64 264 0.498(3) 1.261(17) 1.589(35) 
Asqtad (am,.,d/. = 0.01/0.05), DWF (am„,d = 0.0138) 
28^ X 64 550 0.353(2) 1.191(19) 1.533(27) 
Domain Wall Fermions (DWF) 
m„,d/m, = 0.004/0.03, a'^ = 2.34(3) GeV 



32^ X 64 1452 0.297(5) 1.27(9) 



1.455(17) 



TABLE L Ensembles and parameters used in this work. We 
give in the first column the lattice size, in the second the 
statistics, in the third, fourth and fifth the pion, nucleon and 
A mass in GeV respectively. 



the FFs with small statistical errors. In addition, we 
obtained the FFs using dynamical domain-wall valence 
quarks matched to staggered sea fermions [l^. For the 
computation of the A axial charge we also use Nf = 2-1-1 
domain wall fcrmions poj corresponding to a pion mass 
of about 300 MeV in order to perform the chiral extrap- 
olation. In all cases the u and d quarks are degenerate 
whereas the mass of the strange quark in the dynami- 
cal simulations is set to its physical mass. For the pion 
masses of these simulations the A is stable. 

Lattice results on the A axial form factors: In Fig. [T] 
we show the four axial A form- factors, gi, g^, hi and 
/13, as a function of the momentum transfer Q^. As can 
be seen, results obtained with the mixed action are in 
agreement with quenched results for gi and g^. A similar 
conclusion holds for hi and h^ albeit with much larger 
statistical errors in the case of the mixed action approach 
that wc therefore omit from the plots for clarity. The 
value of the matrix element at = is connected to the 
axial charge defined by (A++|yl3 |A++) - (A-|^3 jA") = 
GaaA^m El- At g2 = this is Gaa = -3gi(0). 

Pseudoscalar FFs and the Goldberger-Treiman Rela- 
tions. The A axial charge enters in baryon xPT ex- 
pressions of many important quantities such as the axial 
charge of the nucleon. Many phenomenological results 
rely on this value, which is usually treated as a fit pa- 
rameter to be determined from fits to experimental or 
lattice data. It can be related to the ttAA coupling via 
the Goldberger-Treiman relation. In Ref. 2l| symmetry 
arguments in a quartet scheme where N^, iV* , and 
A_ form a chiral multiplet, lead to the conclusion that 
7rA±A± couphngs (with like-charged As) are forbidden 
at tree-level. Quark-model arguments [13| suggest that 
the GttAA = (4/5)GTArAr. Clearly a non-perturbative cal- 
culation within lattice QCD of this coupling, as presented 
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FIG. 1: Results for the axial form-factors, gi, gs, hi and /13. 
Noisier mixed-action results are consistent and omitted for 
clarity for the hi and hi FFs. 



in this work, provides valuable input to phenomenology. 

Partial Conservation of the Axial Current (PCAC) 
when applied to the hadronic world leads to valuable 
phenomenological predictions such as the Goldberger- 
Trciman (GT) relation, originally derived for the nuclcon 
state. Similarly, a non-diagonal GT relation, applicable 
to the axial A^ — A transition is formulated and relates the 
axial A^ — A coupling ca to the ttA^A effective coupling. 
PCAC on the hadron level reads: = /.ttTO^tt" . In 

the SU(2) symmetric limit of QCD with denoting the 
up/down mass, the pseudo-scalar density is related to the 
divergence of the axial-vector current through the axial 
Ward-Takahashi identity (AWI) ^^'A1;^ = 2m,P" . Tak- 
ing matrix elements of the LHS of the AWI identity in 
A states we can define two Lorentz-invariant ttAA form 
factors, G'7rAA('?^) and i?7rAA('Z^) factoring out the pion 
pole as dictated by PCAC 



(A(p/,s;)|p3|A(p„s,)) = -2^( 



(4) 



Matrix elements of the AWI identity, (A|a^yl^|A) = 
2TOq(A|P|A) now lead to a matrix equation, satisfied at 
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FIG. 2: The pseudoscalar A form factors CtAA and H-^aa 
for the quenched and dynamical ensembles. 
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where r = —(7^/(2771^)^- We display the pseudoscalar 
FFs in Fig. [5] The results using dynamical quark simu- 
lations have increased statistical errors and are consistent 
with the quenched results. The quark mass mg, extracted 
from the axial Ward-Takahashi identity, and /tt, calcu- 
lated from the pion-to- vacuum amplitude, are taken from 
Ref. [lit . Assuming pion pole dominance, which is con- 
sistent with our lattice results, that indeed 33 and hi are 
of the same order and diverge with the pion pole while h^ 
diverges with (pion pole)^, we obtain a pair of GT-type 
relations, valid at finite q^ , 

/ttG^AaI?^) ^ mAgiiq^) , UH^AAiq'^) ^rriAhiiq'^) .(6) 

The validity of the GT relations is examined 
by evaluating the ratios /7rG7rAA(Q^)/'7iA5i('9^) and 
firHy^AAiQ^) /'mAhi{Q^) as shown in Fig.[3l For the for- 
mer ratio for which statistical errors are smaller the be- 
havior is similar to the one obtained for the pseudoscalar 
nuclcon and A^ — A couplings G-^nn and G-^na [3 for 
the same ensembles, namely for ~ O.SGeV^ the lat- 
tice data show agreement with unity. We expect that the 
behavior at low Q"^ will be affected by pion cloud effects 
as the mass of the pion decreases towards the physical 
point. 

A axial charge and combined chiral fits. Having, for 
the first time, a set of lattice results for the axial nucleon 
charge , the axial TV — A transition coupling, C5 
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FIG. 3: The Goldberger-Treiman ratios from Eq. (O. 



and the A axial charge, allows us to perforin a combined 
fit to all three quantities usin g he avy baryon xPT in the 
small scale expansion [l^ [23I |24| . The combined fit has 
seven free parameters, namely the values of the three ax- 
ial coupling constants of the nucleon, the A'' — A and the 
A, three parameters related to the mj-terms in the chiral 
expansions of gA, and G/\p. and a constant entering 
the chiral expression of C5 [2J] . As can be seen in Fig. SI 
lattice data for all three observables are approximately 
constant within the mass range considered. The best fits 
are shown by the bands that take into account the statis- 
tical errors of the lattice results. As have been observed 
in all recent lattice studies, the physical value of qa is 
underestimated and this combined fit does not provide a 
possible resolution to this puzzle. Having lattice results 
at pion masses below 300 McV will be essential to check 
the validity of these chiral expansions. 

Conclusions. We have presented the first calculation 
of the axial- vector and pseudoscalar form factors of the 
A using lattice QCD. From the most general decompo- 
sition of the axial- vector and pseudoscalar vertex we de- 
rived two Goldberger-Treiman relations whose validity is 
satisfied at the same level of accuracy as that found for 
the nucleon case [l^. At zero momentum transfer the 
A matrix element yields the phenomcnologically impor- 
tant A axial charge, which in this work is obtained for 
pion masses in the range of about 300 McV to 500 McV. 
As in the case of the nucleon axial charge, it shows a 
weak dependence on the pion mass in this mass range. 
Using lattice results for the axial nucleon charge, the ax- 
ial N to A transition coupling and the A axial charge 
we performed, for the first time, a combined fit to all 
three quantities that provides a reasonable description 
to the lattice results. However, these state-of-the-art lat- 
tice results and chiral perturbation calculations, yield a 
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FIG. 4: Combined chiral fit: (a) Nucleon axial charge, qa, 
fitted to lattice data obtained with Nf = 2 twisted mass 
fermions (TMF) The physical value is shown by the 

asterisk; (b) Real part of axial N to A transition coupfing 
C5(0) (c) Real part of A axial charge Gaa = — 35(1(0). 



value for the nucleon axial charge lower than its exper- 
imental value. Such discrepancies between lattice and 
experimental results are seen in several key hadronic ob- 
servables [i^l calling for high accuracy lattice calculations 
with pion mass below 300 MeV in order to to gain insight 
in the chiral behavior of these fundamental quantities. 
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